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ABSTRACT 

We present a matrix solution to the full equations of statistical equilibrium that give 
£3 ' the energy level populations of collisionally-excited ions in photoionised gaseous nebu- 

lae. The rationale for such a calculation is to maintain a parity between improvements 
in the quantum-mechanically evaluated values for collision strengths and transition 
probabilities from the Iron and Opacity Projects on the one hand, and 3D photoion- 
isation codes such as MOCASSIN and astrophysical software for producing nebular 
diagnostics such as the Nebular package for IRAF, on the other. We have taken advan- 
ce \ tage of the fact that mathematics programs such as MATLAB and Mathematica have 
O |. proven to be very adept at symbolic manipulation - providing a route to exact solutions 

for the n-level ion. In particular, we have avoided the substitution of estimated values. 
We provide the matrices for the 5-level ion as an example and show how the equations 
faithfully reduce to the exact solution for the 3-lcvcl ion. Through the forbidden line 
\ ratio i?23, we compare the exact solution with a) that obtained from the observed 

emission of the spherical planetary nebula Abell 39, b) 3D Monte-Carlo photoionisa- 
tion modelling of the same nebula, c) the approximate 5-level program TEMDEN and 
d) the exact 3-level ion. The general solution presented here means that programs for 
the calculation of level populations can obtain solutions for ions with a user-specified 
number of excited levels. The use of a separate and updatable database of atomic and 
ionic constants such as that provided by NIST, means that software of more general 
• application can now be made available; particularly for the study of high excitation 

objects such as active galactic nebulae (AGNs) and supernovae (SNs) where higher 
excited levels become significant. 
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1 INTRODUCTION progress in observational astrophysics has provided empiri- 

cal tools to estimate the values of n e (Menzel, Alter & Hebb 
Essentially 3 parameters fully determine the physical nature ,„.,, ,„ TTrT . , „ , . , rl . „ ,„ r ^ 

, , , , , . . 1941) and T e (Aller, Ufford & Van Vleck 1949; Seaton 1954) 

of an ionised nebula: the electron temperature (T e ), the elec- , . . „. . , , m , 

. , ... and various indicators of ionic abundances (Perez-Montero 

tron density (n e and the level populations of the ions (Ni). „ „, ,. ., , ,, 

, , , . . ,. & Diaz 2005). Observations provide access to the properties 

From the latter, ratios of emission rates and emission line , . . . ,. . T 

of the ionised gas through their emission line spectra. In par- 
mtensities can be calculated. The 3 parameter larmly forms ,. , , . ,. „ „ , . , , i i ii 2 

J ticular, certain ratios of forbidden lines emmited by the p 

a closed set whereby knowledge of any two of them allows .3 r . . , . ,, , , 

and p conhguration 10ns present m the nebular gas, allow 
for a determination of the third. In principle, all other sec- , , . , , n n , 11 / 

the integrated electron density and electron temperature (re- 
ondary physical quantities such as lomsation parameters, ,. , n , , , , . , . , . 1 , 

, , , rr 11 spectively) to be estimated point to point across projected 

ionic abundances and effective temperatures can then be . „ rr , . , , ,. .. . , , n . r 

. , , TT , images, furthermore, sums ol torbidden lines emitted by dif- 

calculated. Until now the exact functional relation between „ , r . . , . , , , • i- , \ 

_ , . terent stages ol lomsation ol the gas (abundance indicators) 

the parameters has evaded astronomers. In this paper, we , , r , . , , , . .., , „. ., TT 

. have been lound to weakly correlate with metalhcity. How- 

present a solution to this problem. The context ol great , , . r 1 1 1 

ever, the underlying reason lor such correlations has yet to 

be totally resolved. In this paper we show that, since the 

emission line intensities depend on the level populations of 
* michael@damir.iem.csic.es . . . . . , , , . , . . 

, . ^. the 10ns which m turn depend on the local electron density 

j jvm@iaa.es 1 ■' 
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and temperature, then such correlations simply reflect an 
empirical relation between T e , n e and Ni of the form, 

f(n e ,T e ,Ni) = 0. (1) 

Using linear algebra, we have re-cast the equations of 
statistical equilibrium for the level populations of the 
ions in the form of a well-posed matrix equation that we 
subsequently solved using symbolic manipulation to obtain 
the exact form of the function /. 

The early work in this field was performed by Men- 
zel(1941, 1962), Aller(1949), Seaton(1954, 1975) and 
Osterbrock(1967) who obtained initial estimates of the level 
populations of 3-level ions. The full solution for the 3-level 
ion was finally worked out by Seaton(1975) and a first order 
approximation the 5-level ion is currently incorporated 
into astrophysical software such as TEMDEN (De Robertis 
et al. 1987). We use these results as well as detailed 
spectrophotometric observations and 3D Monte Carlo 
photoionisation modelling of the spherical planetary nebula 
Abell 39 to validate and compare with the exact solution 
we have obtained. The theoretical context for the physics of 
collisionally-excited ions is the thermal equilibrium present 
in ionised gases. In the next section we highlight the role 
played by collisionally-excited emission lines (and therefore 
of the n-level ion) to the thermal balance of ionised gaseous 
nebulae. 



2 THERMAL EQUILIBRIUM 

The temperature in a static nebulae is fixed by the equi- 
librium between heating by photoionisation and cooling 
by recombination, free-free radiation (bremsstrahlung) and 
line radiation. In ionisation equilibrium, photoionisations 
are balanced by an equal number of recombinations of 
thermal electrons. The difference between the mean energy 
of a newly created photoelectron and the recombining 
thermal electron represents the net gain in energy of the 
electron gas per ionisation process. In equilibrium, this 
net gain is balanced by the energy lost by bremsstrahlung 
and radiation produced by collisional excitation of bound 
levels of abundant ions that subsequently emit photons 
leaving the nebula. For radiatively-ionised, low density 
plasmas, the relaxation times of relevant physical processes 
are such that we can usually regard a given volume element 
as being in a steady state. If we can neglect mechanical 
effects, particularly energy lost or gained by expansion 
and compression, we can impose the condition of radiative 
equilibrium (or energy balance) such that the energy 
absorbed in each volume element then equals the amount 
emitted (Spitzer 1948; Osterbrock & Ferland 2006). 

Menzel and co-workers first set up the equation of en- 
ergy conservation (Menzel et al. 1941) whereby the energy 
absorbed by photoionisation of the gas (G) is balanced by 
the energy liberated in capture and subsequent recombina- 
tion events (Lr), free-free Bremsstrahlung emission (£//) 
and energy emitted in collisionally-excited radiative cooling 
(La), 

G = L R + L ff + L c . (2) 



Closed forms for G, Lr and L// can be found in standard 
texts on the physics of ionised gaseous nebula (e.g. Aller 
1984; Osterbrock & Ferland 2006). What interests us here 
in the present work is the form of Lc, the cooling due to 
the emission of collisionally-excited lines (CELs) and their 
predominance over the estimation of ionic tempertures. At 
densities higher than the critical density N c , 

j<i j^i 

collisional de-excitation becomes important and the cool- 
ing rate at a given temperature is decreased. Hence, for 
n e < N c (i), collisional de-excitation of level i is negligible 
- the thermal regime of the vast majority of planetary 
nebulae. 



3 ENERGY LOSS BY CELS 

In many gaseous nebulae, mechanical or magnetic energy 
from hydromagnetic waves is dissipated in the gas and 
there can exist regions where the mean temperature is 
raised such that atomic levels are excited through collisions 
with electrons. Furthermore, in the presence of sources of 
high energy UV photons such as those produced in stellar 
atmospheres or in the radiation fields of stellar clouds, 
thermalised electrons are capable of ionising the low-lying 
levels of the ions of heavy elements. As these ions return to 
their ground states, they emit permitted line photons (that 
fall in the UV domain observable with satellite telescopes). 
However, atomic term structure, elemental abundances 
and ionic concentrations are such that in gaseous nebulae, 
very few collisionally excited permitted lines are observed 
in optical spectral regions, particularly in the low density 
plasmas of many nebulae and HII regions. The most 
frequently observed optical lines come from ions excited by 
impacts with electrons - the forbidden lines that violate the 
Laporte Parity Rule (Aller 1984). For these lines, we need 
to know the transition probabilities (Atj) and the collision 
strengths (flij) for the excitation of metastable levels. 
Forbidden lines cannot be observed in the laboratory and 
so we have to rely entirely on theory to determine Aij and 
Oy. Observational checks on certain line intensity ratios are 
sometimes possible, but for most of the transitions on which 
we depend for nebular plasma diagnostics, observational or 
experimental checks are not available. In the next section 
therefore, we rigorously solve for the populations of the 
n-level ion that give rise to emitted forbidden lines. 



3.1 The n-level ion 

Assuming that level populations Ni are determined by spon- 
taneous emission (with transition probability Ay per sec- 
ond), i.e. we note that the only possible radiative transitions 
are downward ones resulting from impacts with thermal elec- 
trons, then the formal general equations of statisitcal equi- 
librium for an ion with excited levels (i) above ground having 
excitation rate coefficients (qji) and de-excitation rate coeffi- 
cients (qij) to the other levels j ^ i are given by Aller(1984) 
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equation 5-32), 



with 



Ei are the 



Nin e qij + NjAjj 

j<i 



(1) 



jjii j>i j^ti 

subject to the condition that the sum over all level popula- 
tions (i.e. over all stages of ionisation) is equal to the total 
number density (TV) of the ions [cm -3 ], 



E 



Ni = TV. 



(5) 



Here we have adopted the excitation potential convention 
Ei > Ej such that the transition i — » j corresponds to 
de-excitation. The balance is between terms that contribute 
to the population of an ionised level and those that subtract 
from it. The first term on the left hand side therefore 
represents collisional de(excitations) to i from other levels 
(j 7^ i) while the second term includes radiative transitions 
to i from upper levels (j > i). The left hand side therefore 
depends on the populations of levels other than i. The 
first term on the right hand side represents collisional 
de-excitation of other levels by electrons from i while the 
second term represents downward radiative transitions from 
level i to lower levels (j < i). The right hand side therefore 
depends on the population of level i. 

The next step is to note that qij is related to the 
collisional excitation rate dj [s~ ], 

js- 

Cij — Qij = l\ X 

UJi 



(6) 



with K = 8.629 x 10" 6 , Seaton variable x = n e /T e 1/2 and 
where flji is the velocity-averaged collision strength (Seaton 
1968), 



Wme--'™*(JLr), (7) 



for colliding electrons having inital kinetic energy 
E = mu 2 /2. Here, T e is the electron temperature [K] 
and n e is the electron density [cm -3 ]. The collision 
strengths must be calculated quantum-mechanically and 
consist, in general, of a part that varies slowly with energy 
superimposed with resonance contributions that vary 
rapidly. Early attempts at calculating them using the 
Born-Oppenheimer approximation turned out to violate 
the Mott-Bohr-Peierls-Placzek conservation theorem (Aller 
1984), and it wasn't until the late 1950s and early 1960s 
that Seaton successfully calculated reliable cross-sections 
using the "Exact Resonance Method" and the "Distorted 
Wave Method" (Seaton 1968, 1975). With the development 
of the close-coupling approximation, consistency arguments 
and, most importantly, the fact that integration over a 
broad Maxwellian distribution of electron energies tends to 
dramatically smooth out variations, Qji is now known to 
be fairly insensitive to temperature (Osterbrock & Ferland 
2006). 



excitation potentials of the levels. We see that, through the 
collisional de(excitation) rates and their dependence on the 
Seaton variable x(n e ,T e ) and €ij(T e ), the level populations 
(and therefore all quantities derived from them such as 
emission line ratios) have the dependency on both electron 
density and electron temperature given in equation 1. In 
contrast, the radiative transition probabilities Ai t j are 
independent constants being inversely proportional to the 
occupancy lifetimes of the upper level. We note in passing 
that for many nebulae the Seaton variable « 10 2 . 

Replacing all n e qji and n B qij terms in the equations 
of thermal equilibrium for the excited levels by the expres- 
sions above for the collisional excitation and de-excitation 
rates and dividing through by Kx / we obtain, 

j<i 



Kx 



^m£i= . (io) 



We therefore have a set of simultaneous equations in n un- 
knowns for the (n — 1) excited levels above ground (i — 1), 
supplemented by the total ion density condition Ni + N2 + 
■ ■ ■ + N n — TV. Taken together as a set, they have the math- 
ematical form, 



(11) 
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1 


1 


1 




TVi 




TV 


«21 


CK22 


«23 


0?2n 




TV 2 







«31 


«32 


Q33 


Q3„ 




TV 3 







"41 


C«42 


«43 


C<4n 




TV 4 







C<51 


CS52 


«53 


C<5n 




TV 5 







dnl 


a n 2 


Ctn3 














where the a%j are the coefficients of the TVi. The first row 
with OL\j — 1 reflects the total ion density condition (equa- 
tion 5). Introducing the matrix A of coefficients ajj, the 
vector y of level populations TVi and the vector b for the 
right hand side, we can write the simultaneous equations as 
the linear matrix system, 

Ay = b. (12) 
Provided that | Ji| 7^ then we can solve for y, 

y = A- 1 h (13) 

using symbolic manipulation software such as MATLAB or 
Mathematica. The sought-after level populations are identi- 
cally given by, 



TVi 


= Vi 


= (A- 


x )nN 


TV 2 


= J/2 


= (A- 


%iN 


TV 3 


= J/3 


= (A- 


%iN 


TV 4 


= J/4 


= (A- 


%iN 


TVs 


= J/5 


= (A- 





The collisional de-excitation rate, 

Cij = n e qji, (8) 
is related to the collisional excitation rate through, 

Cji = Cij €ij, (9) 



= y n = (A 1 ) n iN. 



(14) 



Finally, the collisionally-excited radiative cooling rate (Lc), 
central to the overall thermal balance of gaseous nebulae, 
can then be calculated from the level populations and is 
given by, 
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(15) 



where hfij is the difference in energy level potentials. 



Equations 10-14 allow the astrophysicist to model (exactly) 
an ion ol as many levels as required, provided |A| 7^ such 
that the matrix inverse A -1 exists. In the next section we 
make this more explicit by providing the matrix equations 
for the 5-level ion. 



3.2 The 5-level ion 

Ions having p 2 , p 3 and p 4 electron configurations all have 
5 low-lying energy levels. For such ions, collisional and 
radiative transitions can occur between any of the levels 
and excitation and de-excitation cross-sections as well 
as collision strengths exist between all pairs of levels. A 
central assumption that has historically been made up until 
now is that only these 5 levels are physically relevant to a 
calculation of the observed emission lines of the ion. The 
justification (Osterbrock & Ferland 2006) is that higher 
levels in these ions are not significantly populated through 
collisions, recombinations or other mechanisms. However, 
for strongly ionised nebulae in AGNs or SNs for example, 
the higher levels may indeed be significant to the emission 
line spectrum. 

For a 5-level atom in a steady state with E$ > E 4 > 
E3 > E2 > Ei, the total number density condition for each 
ion species is given by, 



N 1 +N 2 +N 3 +N 4 +N 5 = N. 



(16) 



For the four excited levels above ground, the equations of 
statistical equilibrium (equation 10) with statistical weights 
LOk = (2fc + 1) give rise to the following four exact level 
population equations: 

f^2l£21 



iVi- 



\ Kx 



+ D. 



+ N 4 
+N 5 



',32 + ^42 


+ ^52 


5 




(A* 


^42 £24 


\Kx 


9 


(*» + 
\Kx 


^52^25 


11 



= 



(17) 



-Na 



f2l3£31 , , T ^23^32 
iVl h 1\2 = 

3 5 

A31 + ^32 , «13 + ^23 + ^43 + ^53 



Kx 7 



f2l4£41 , Ar ^24£42 . », ^34£43 
iVl h iV2 h JV3 - 



(18) 



-N 4 



5 " " 7 
'An + A 42 + A 43 



-N 4 



Kx 

' S2l4 + ^24 + ^34 + f254 



(19) 



Ni 



^15^5 



^45 £54 



5 • - " 7 '9 
Asi + A52 + A 53 + A 54 



-N 5 



Kx 

' Ol5 + 5^25 + ^35 + "45 

v n 



(20) 



The elements ay of A are then: 
an = ai2 = Q13 = Q!14 = ai5 



«21 = 

Q22 = 

Q23 = 

Q24 = 

«25 = 

Q31 = 

Q32 = 

«33 = 

«34 = 

C*35 = 

Q41 = 

Q42 = 

«43 = 

Q44 = 

«45 = 

«51 = 

Q52 = 

C*53 = 

«54 = 

C*55 = 



^12621 

3 

A21 
Kx ' 



fil2 + ^32 + ^42 + f2 5 



Kx 


+ 


7 


A 4 2 




^42 £24 


Kx 


+ 


9 


A 52 


+ 


^52 £25 



fil3£31 
3 

^23^32 



(A31 + A3 



+ 



Kx 

f2l3 + Q23 + ^43 + ^53 
7 

A 4 3 £l 4 3 £34 

~Kx~ 9 

^4-53 ^53 £35 

Kx 11 



fil4£41 

3 

24 e 4 2 
5 

^34£43 
7 



+ 



(Am + A 42 + A 43 ) 
Kx 

Ql4 + f^24 + ^34 + f^54 

9 

Ar, 4 fls4£45 

~Kx~ 11 



^15£51 



3 

^25 £52 

5 

^35£53 

7 

$145654 
9 



(A51 + A52 + A 53 + A 54 ) 
Kx 
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il 



(21) 



It is at this point that we can identify the function / in 
equation 1. The population of the jth level Nj is given by, 

Nj = (A'^iN. (22) 

Using the notion of the matrix of cofactors Cji known as an 
adjugate matrix such that, 

111 



(23) 



then the level populations can be written in the form, 
C , J i^(r e ),x(n e ,T e ))iV 



N 



or alternatively, 

iV i =/(n e ,T.) 



(24) 



(25) 



as implied by equation 1. We have deliberately not inserted 
the values of Aji, Qji and up until now as they are 
estimated values that will evolve with the accuracy of quan- 
tum mechanical calculations currently being performed by 
the Iron and Opacity Projects. In its present form, the 
matrix equation Ay — b with coefficients ay is exact. This 
is the closed form for the set of simultaneous equations 
for the 5-level ion. We used the symbolic manipulation 
toolbox of MATLAB to find the inverse matrix A' 1 and 
the level populations. The full solution includes terms 
up to (Kx) 4 and is too big to reproduce here. In the 
online supplementary material for this paper, we provide 
a MATLAB m-file whose results can be included into 
astrophysical codes. In general, an n-level ion will contain 
terms up to (Kx)"" 1 . The solution for the 3-level ion in- 
cludes terms only up to (Kx) 2 and is much more compact. 
Furthermore, its exact solution is known (Seaton 1975). 
This then provides a vital algebraic check on our solu- 
tion to the 5-level ion as we demonstrate in the next section. 



3.3 The 3-level ion 

Supressing all terms with indices equal to 4 and 5 reduces 
the general solution for the 5-level ion to that of the 3-level 
ion. In doing so, for ions of type p 2 such as [OIII] as shown in 
figure 1, this amounts to removing the fine splitting physics 
of the 3 P level. As a consequence, the excitation potentials 
E2 and E3 of the 3-level ion correspond identically with E<± 
and E5 of the 5-level ion while Ei is taken to be ground. 

The exact level populations for the 3-level ion are then given 
by: 



3iV 



35^21^31 + 35^21^32 
+Kx^5A 2 lQl3 + 5^21^23 + 7^31^12 
+7A 31 Q 32 + 7^32^12 + 7^32^32 - 7^32^23^32 
+K 2 x 2 ( ^12^13 + ^12^23 + ^13^3 



+^23^32 — £723 632 5732 623 ) (26) 



4363.2 



D 



233 1 .4 



5 E 3 = 5.354 eV 



5006.8 | 
3 I 4958.9 



493 1 .2 



2321,0 



4 E,= 2.514eV 



-3 E 3 = 0.038 eV 
* — 2 E 2 = 0.014 eV 
1 E,= eV 



Figure 1. The 5-level [OIII] ion. Radiative transitions marked 
by arrows have their central emission wavelength in A 



N 2 



- 7jfic ^4.31 fil2 621 +^32^12^21 + A^Vi^il 

+K 2 X 2 ^fii2C21 fil3 + 57i262l 57 2 3 + 57i 3 63i 57 3 2623 



N 3 = 



7N 
~D 



Kx\ 5^21^13631 



+K 2 X 2 ^12621^23632 + 57l 257 l363l + 57l 3573 l 5732 



(27) 



(28) 



with the denominator D term, 

D = 105^21^31 + 105^21^32 
+KX 15A 2 lfil3 + 15^21^23 + 21^31^12 + 21^31^32 
+21^32^12 + 21^32^32 - 21^32^23632 + 35^21^13631 
+35^31^12621 + 35^32^12621 + 35^32^13631 

+K 2 X 2 357i 2 57i 3 + 3^12^23 + 357i 3 57 3 2 

+3^23^32 — 35723632^32 623 
+5012621^13 + 5^12621^23 + 5^13631^32623 

+75712^13631 + 75112621^23632 + 7^13631^32 . (29) 

3.4 Line ratios of the 3-level ion 

The intensity of a radiatively-emitted line produced by col- 
lisional excitation of ions in a nebula can be written down 
directly from the level populations by noting that for any 
emission line, the emission rate j(i, k) of line photons re- 
sulting from a downward transition i — > k is given by, 



Anj(i, k) = A ik N(X l )Nihv ik 



(30) 



where X 1 is the ion species e.g. 2+ and hvik corresponds 
to the energy difference AE between excitation potentials 
Ei and Ek- The number of ions relative to ionised Hydrogen 
is given by, 
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N(X') I(i,k)j(H p ) 



(31) 



N(H+) I(Hp)j(i,ky 

As stated in the introduction, there are standard methods 
that use optical emission line intensity ratios to measure 
electron temperature (Menzel et al. 1941) and electron den- 
sity (Aller et al. 1949; Seaton 1954), including the simulta- 
neous determination of both temperature and density (Din- 
erstein, Lester & Werner 1985). Comparisons have also been 
made between various density indicators (Copetti & Writzl 
2002) and temperature indicators (Taylor & Vflchez 2007d) 
as well as methodologies for the comparison between the- 
ory and observation (Seaton 1960; Taylor & Diaz 2007a). In 
particular, the line ratio R23, 

7(A4959) + /(A5007) 



R23 



(32) 



/(A4363) 

is a well known temperature diagnostic in ionised gaseous 
nebulae (Osterbrock & Ferland 2006) and is based on the ion 
[OIII] (Pagel et al. 1979) whose first 5 levels above ground 
are shown in figure 1. In taking the ratio of emission lines, 
the number densities of the ion and ionised Hydrogen cancel 
leaving an expression in terms of atomic and ionic constants 
only. So, for example, for the 3-level ion, R23 is given by, 

j(2, 1)5007 +j(2, 1)4959 



-R23 



3(3,2)4 



N 2 AE 21 
N3AE32 



A21 (5007) + ^21 (4959) 



^32 (4363) (33) 

Many such line ratios of the ions of He, C, N, O, Ar, S, Ni 
and CI are used as diagnostics of the physical conditions in 
photo-ionised gaseous nebulae (application of the formulae 
derived here for the n-level ion to a large collection of such 
diagnostics is in preparation(Taylor & Vflchez 2007b). 

For the purposes of the present discussion, we con- 
sider the ratio of the level populations N3 and N2 (although 
nowadays the inverse N2/N3 for the ratio of nebular to 
auroral lines is customary) so as to compare our exact 
solution for the 3-level ion with the solution obtained by 
Seaton(1975): 



IN 
D 



Kx[ 5^4 2 if2i3£3i 



+K 2 X 2 I r2l2£21^23£32 + fil 2 f2l 3 e31 + ^13^31^3 



N3 

N 2 



D 



7Kx[ A 31 Q 12 €21 + ^4-32^12621 + ^32^13631 

+K 2 X 2 ^Qi2<^\U,: i + ih^-jiih.: + 12 1.-.<-. i S. »:,.:< _m 

Extracting the term (^13631) from the numerator and 
(f2i2£2i) from the denominator and calcelling like terms we 
get the result, 

iV 3 _ gig 631 A 2 l c [l+^3] 
N 2 fii2 £21 A32 [l+*3Ltf 2 ] 

where, 

1 

C = 



* 2 = 



1 + 

c( 



+ 



^13 g 31 



^32 ^12 e 21 
^12621^13 + f2l2£2lfi23 + ^13e31^32£23 ' 
7^12621 



, T , fil2£2lfi23£32 + ^12^13£31 + f2i 3 e 3 if232£23 ,„_ x 

*3 = ^— . (35) 

Oi2l3e31 

This is the same result as that obtained for the exact 3-level 
ion (Seaton 1975: equation 1.11). 



4 THEORETICAL, OBSERVATIONAL AND 
EXPERIMENTAL CHECKS 

As with all new theoretical results, it is vital to perform 
checks. In the previous section we have verified algebraically 
that the equations for the 5-level ion (calculated from the 
full n-level ion equations), when reduced to the case of 
the 3-level ion, reproduce earlier results. In this section, 
we perform two additional checks. By considering the 
dimensions of the physical variables, we perform a check 
on the dimensional homogeneity of the solutions. Then, 
by inserting the most up-to-date quantum-mechanically 
calculated values for Aj%, tij and Qji, we compare the 
values of the forbidden line ratio R23 from the exact 
equations for the 5-level ion with those obtained from: a) 
observations of the benchmark planetary nebulae Abell 
39, b) 3D Monte-Carlo photoionisation modelling of A39, 
c) numerical approximations provided by the nebular 
diagnostic software TEMDEN (De Roberties et al. 1987), 
and d) the exact 3-level ion. 



4.1 Dimensional consistency of the equations 

There are 7 units in the SI: mass M[Kg\, length L[m], time 
T[s], temperature 0[K], electrical current I [A], concentra- 
tion N[mol] and light intensity J[cd]. In all equations of 
theoretical and empirical science, dimensional homogeneity 
for any set of physical variables {xi} must satisfy, 



= 1. 



(36) 



for causality to be satisfied (Taylor et al. 2007c). The di- 
mensions of the various physical quantities and parameters 
presented in our derivation are, 



such that, 



K] = [Ni] 
[qji] = [qij] 
[n e Niqji] = 

[Aji] 
[x(n e ,T e )] 

[K] 

= [fyi] = [ay] 



j3 T -\ 

L -3 T -1 

T -l 

L -3 e -l/2 

L 3 r _1 e 1/2 
1 



Aji 




' Kx 


IKx. 




Aji ^ 



= [C] = [9 2] = [* 3 ] = 1 



All equations were thoroughly checked for dimensional 
homogeneity with these dimensional relations and we did 
not find any inconsistencies. 



© 2006 RAS, MNRAS 000, [Lj{8] 



Exact solutions for the populations of the n-level ion 7 



4.2 The forbidden line ratio R23 



5 CONCLUSIONS 



We have used observations of the spherically-symmetrical 
planetary nebula Abell 39 to make comparisons with our 
theoretical predictions. Although discovered by George 
Abell in his 1957 survey of the Southern Hemisphere (Abell 
1966), A39 has only recently been subject to thorough spec- 
trophotometry(Jacoby, Ferland & Korista 2000) and 3D 
photoionisation modelling(Taylor et al. 2007e). From the de- 
tailed measurements made by George Jacoby and co-workers 
at Kitt Peak, the redenning-corrected line ratio R23 averaged 
over the whole nebula, given by equation 28 equals, 



R23 = 



±0.182 + 11.31 ±0.524 
0.24 ±0.019 



: 64.24 ± 6.58 



(37) 



with all fluxes I(X) relative to I (Hp). Using the 3D Monte- 
Carlo photoionisation code MOCASSIN (Ercolano et al. 
2003) and a best-fit model to the observed spectrum(Taylor 
et al. 2007e), the following value, well within the observa- 
tional error, was obtained, 



We have shown that it is fairly straight-forward to gener- 
alise the calculation of ion level populations from the formal 
general statistical equilibrium equations. The matrix formu- 
lation presented here together with the algebraic power of- 
fered by symbolic manipulation software means that it has 
been possible to obtain an exact solution for the n-level ion. 
As a check on this result we have shown that, when indices 
for levels 4 and 5 for the 5-level ion are supressed, the exact 
solution for the 3-level ion is recovered. Iterations to higher 
level ions are straight-forward and we hope that these results 
will make the incorporation of more ionic transitions into ex- 
isting astrophysics software easier for the study of very high 
ionisation objects such as AGN and SN. The 3-fold compar- 
ison between observations, theory and simulations appear to 
be consistent. In Paper II we will investigate a set of line ra- 
tio diagnostics for low, medium and high zones of ionisation 
based on the exact solution obtained here. 



3.51 + 10.47 
= 021 = 



(38) 



In Table 1, we calculate the value of R23 for the 5-level 



D N 4 r AS43A43 + AE 4 2Ai2 

R23 = w 5 



AE 54 A 5 4 



(39) 



using our exact theory, inserting the most up-to-date 
values of the atomic and ionic constants Aij, Qi t j and tji 
(Osterbrock & Ferland 2006 and references therein), for 
representative values of electron density and temperature 
deduced empirically from the observations and simulations 
of A39. In addition, we provide results at a lower electron 
temperature of 10,000i(" as a more extreme test case. The 
"theoretical" value of R23 (calculated from the exact the- 
ory) is compared with that obtained from a) observations, 
b) 3D simulations using MOCASSIN, c) the first order 
5-level ion of TEMDEN and d) the exact 3-level ion. 

At electron densities of n e = 30cm -3 all of the re- 
sults are within the observational standard error apart 
from the low temperature case (T e = 10000) where, even 
here, there is consistency between the exact 5-level theory, 
TEMDEN and the exact 3-level theory. In the case of 
the best-fit 3D model with MOCASSIN and using a 
3-shell density profile having an average electron density 
of < rt e >= 11.94cm -3 , the results are more dispersed 
with the 3D code producing a value of R23 closer to the 
observations. As we have mentioned, this may be due to 
uncertainties in the values of the atomic constants Aij, 
and eji. What is interesting is that the value of R23 
from the exact 5-level ion treatment is consistently higher 
than both the exact 3-level ion and also the first order 
approximation to the 5-level ion of TEMDEN. The effect 
of including all of the atomic physics provided by the exact 
theory does not appear to play a very significant role at 
these lower excitation stages nor at moderate values of 
plasma density and temperature. 
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Table 1. Observed, simulated and theoretical values of the [OIII] line ratio R23 



R23 

CASE STUDY n e T e Observed MOCASSIN TEMDEN 5-level 3-level 
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(Jacoby et al. 2000) 
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(3-shell density profile) 
(uniform density) 
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n e = 30 



(T e ) = 15896 
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69.91 
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61.02 
69.22 
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61.27 62.37 60.34 
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